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'. Abstract 

We confront tracker field quintessence with observational data. The potentials considered in this 
■ paper include V{(p) oc </> exp(Mp/(/)), exp(Mp/(/)) — 1, exp(/?Mp/(/)) and exp(7Mp/(/)) — 1; while 

'sj- . the data come from the latest SN la, CMB and BAO observations. Stringent parameter constraints 
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are obtained. In comparison with the cosmological constant via information criteria, it is found 
that models with potentials exp(Mp/(/)), exp{Mp/(p) — 1 and exp{'yMp/(p) — 1 are not supported by 
the current data. 
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I. INTORDUCTION 



Research on the cosmic accelerating expansion from both observational and theoretical 
sides has been proceeding intensely for the past dozen years (see jl, 2| for reviews). Within 
the framework of general relativity, the cosmic acceleration requires that the universe's cur- 
rent energy budget is dominated by a energy source, named dark energy, providing significant 
negative pressure density p < — p/3, where p is its energy density. Combined constraints 
from different types of current observation 3|- 6| render the best-fit dark energy equation 
of state, w = p/ p ^ —1 with a few percent uncertainty, assuming w is constant. When 
dynamical w is considered, its current value is still close to —1, with about 10% uncertainty. 
This indicates that the cosmological constant remains the simplest valid realization of dark 
energy but there is still room for other possibilities. 

Phenomenological studies on dark energy can in principle be categorized into two ap- 
proaches. One is to reconstruct general properties of dark energy, the other is to constrain 
models on an individual basis. In the first approach, the evolution of equation of state w{z), 
where z is the redshift, for instance, can be reconstructed using either piecewise parametriza- 
tion, continuous parametrization, or principal component analysis (see 4], j?] and js] for 
recent examples). The joint evolution of w{z) and its time-derivative in units of the Hubble 
time, w' = dw/dlna, can also be reconstructed in comparison with theoretical boundaries 
for various dark energy classes These studies render us general features of dark energy, 
yet the results may vary depending on the parametrization in use. Furthermore, the recon- 
struction technique has been developed to provide diagnostics [lO| and consistency tests for 



dark energy models 13|. While these tests can in principle be taken to falsify models, 
we should be aware of the possible bias carried by the the chosen parametrization 

Orthogonal but complimentary to the reconstruction is the model-based approach, in 
which we directly constrain the parameter space of a dark energy model using observa- 
tional data (see [l^ for the case of pNGB quintessence, for example). As the precision 
of observation advances, this approach is becoming effective. Whereas robust and stringent 
constraints on model parameters can be obtained, a model's validity and its relative merit to 
other contending models can be evaluated via the goodness of fit (GoF) and the information 
criteria. 

In this paper, we take the model-based approach to confront the tracker field quintessence 
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model [15|, [iGj with observational data. In the quintessence scenario, the late time cosmic 
acceleration is driven by a scaler field witch slowly rolls down its potential. As a class of 
quintessence model, the tracker field has an attractive feature that there exists a common 
solution to the equation of motion, extremely insensitive to initial conditions. This feature 
can address the cosmic coincidence problem, that is, the ratio of the dark energy density 
to the matter density must be set to a specific, infinitesimal value in the early universe 
in order to be of the order of one today. The data we use come from the latest Type la 
supernova (SN la) compilation, the cosmic microwave background (CMB) observation, and 
the observation of the baryon acoustic oscillations (BAO). Besides constraining the model 
parameters, we assess the GoF and the model strength in comparison with the cosmological 
constant model. 

II. TRACKER FIELD QUINTESSENCE 
A. Quintessence formalism 



In the quintessence scenario [17j , the late time cosmic acceleration is driven by a dynami- 
cal scalar field slowly evolving in the potential V{(f)). In a fiat Friedmann- Robertson- Walker 
(FRW) universe, the evolution of the scalar field is governed by its equation of motion 

^ + 3i/^+| = 0, (1) 

where the dots denote derivatives with respect to time, H is the Hubble expansion rate a/ a 
(a is the scale factor) given by the Friedmann equation 



Hi 



+ zf + + zf + fi^exp (^3 ^ [1 + w^{z')] 



(2) 



where Pr{z) is the radiation energy density, Pm{z) is the matter energy density, P(f){z) is the 
scalar field energy density, Hq is the Hubble constant and w,f){z) is the equation of state of 
the scalar field. The total fractional energy density today Qtotai = + + ^0 is equal to 
one in a fiat universe. The energy density and pressure density of the scalar field are 

P0 = ^0' + n0), (3) 
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P^ = ^^'-V{<P). (4) 

The equation of state W(j){z) = p^/p^ changes with time and becomes negative when the 
potential is dominant. In the hmit when 0^ ^ ^(0)? the scalar filed has w^{z) ~ —1. 



B. Tracker fields 

Tracker fields are a class of quintessence that address the coincidence problem. In these 
models, a wide range of initial conditions in the early universe evolve toward a common 
solution, called tracking solution, giving the same late time evolution of 0, a nd P d,, and 



allowing the scalar field to drive the cosmic acceleration. The central theorem in [16| states 
that tracking behavior with w^f, < wb, where wb is the equation of state of the dominant 
background component, occurs for any potential in which F = V"V/{V'y > 1 (the primes 
denote derivatives with respect to (p) and is nearly constant |r'/r(y/V")| ^ 1 over the 
range of plausible initial conditions. The range of plausible initial conditions extends from 
Ptf, equal to the initial radiation energy density pr in the early universe down to P(j, equal 
to the current matter energy density pm- This constraint is necessary for the scalar field to 
converge to the tracking solution before the present time. The feature Wcf, < wb means that 
P0 decreases more slowly then the background energy density. Eventually, at late time the 
scalar field density overtakes the matter density and becomes the dominant component. 

A class of potentials that satisfies the conditions that F > 1 and F is nearly constant 
includes the inverse power-law (IPL) potential oc and combinations of IPL terms, 
for example V{(f)) oc exp(Mp/0), where Mp is the Planck mass. Some of these potentials are 
motivated b y p article physics models with dynamical symmetry breaking or nonperturbative 
effects [isll 1241 . The illustration in Fig. [1] exemplifies the tracking behavior and the late 
time tracking solution of these models. Details of the tracker field property can be found 
in 



15|, [l6|, l25|] . Analytical solution to the IPL model has been studied in |26l. |27|. 
In this paper, we analyze the models exemplified in [igI], ^(0) oc 0~", ^(0) oc exp(Mp/0), 
and V^(0) oc exp(Mp/0) — 1. We also analyze the generalization of the last two, V"(0) oc 
exp(/3Mp/0) and V^(0) oc exp(7Mp/0) — 1. All a, (3 and 7 are positive constant. Both 
1^(0) oc exp(Mp/0) and 1^(0) oc exp(Mp/0) — 1 are distinct from the cosmological constant 
without extra parameters being introduced. V(0) oc 0~" and V^(0) oc exp(/3Mp/0) behave 
like a cosmological constant as a and (3 approach zero, respectively. 1^(0) oc exp(7Mp/0) — 1 



does not have a limit as the cosmological constant. 



III. DATA 

We use observational data from SN la, CMB and BAO as described below. 
A. Type la supernovae 

We use the latest SNe la dataset, Union2 compilations, released by Supernova Cosmology 
Project which contains 557 SNe la in the the redshift range 0.02 < 2; < 1.5 [4]. This 



compilation includes supernova data from 



28|-[36|. The dataset provides distance modulus 



which contains information of luminosity distance that can be used to constrain dark energy. 
The distance modulus is defined as following: 

fith{z) = 5/05(10 f ) + 25 = 5/05(10 {Dl{z)) + fio, (5) 



Mpc 

where Dl{z) = HqcIlIz) is Hubble- free luminosity distance. We marginalize XsNia o"^^^ 
the nuisance parameter /io by minimizing it with respect to fiQ. The marginalized XsNia 



is |37|-|39| 



where 

^ = ^ (5/0510 [DL{zi,par)] - Hobs{zi)) C^^ {5logio [DL{zj,par)] - Hobs{zj)) , (7) 
B = ^ (5/0510 [DL{zi,par)] - ^iobs{zi)) (8) 



B. Cosmic microwave background 

The seven-year WMAP results provide "distance prior" that can be used to constrain 
dark energy [3]. Distance prior includes CMB shift parameter R = 1.725 ± 0.018 given by 



R= jQ^Hi{l + z,)DA{z,), (10) 
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FIG. 1: An illustration of the tracking behavior and late time evolution of u'^(z) of the tracking 
solution. The models include V{(j)) oc (/>^", cxp{/3Mp/(p) and exp{'yMp/(p) — 1. is set to be 0.27. 
Ptf, starting from a wide dynamical range in the early universe all converge to the same tracking 
solution at late time. For all the three models, W(j,{z) decreases in time at late time, the smaller 
the model parameter is the more negative value w^{z) can reach. 
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and "acoustic scale" Ia = 302.09 ± 0.76 given by 

, ttDa (-2*) / s 

Ts (-2*) 

where is the redshift of decouphng, Da is the angular diameter distance, and is the 



comoving sound horizon. We use the fitting formula proposed by Hu and Sugiyama [40|: 



1048 [1 + 0.00124(n6/i2)-o.738j ^ g^{Q^hy^] , (12) 
0.0783(^],/^^)-°^^3« 

l + 39.5(l]b/l2)0-763' ^ 

^ 0.560 

1 + 21.1(^^/12)1.81- ^ -I 



The comoving sound horizon is 



r.Jz) = — / (15) 

73 io a2i7(a)v/rT(3rV4^M^ 



where Qb is baryon density and fi^ is photon density. 

We construct Xqmb — ~ ~ ^f''^), where C^j^ is the inverse covariance 



matrix given in 



and Xi = {I A, R, z^). 



C. Baryon acoustic oscillations 

We use BAO data from the joint analysis of Two Degree Field Galaxy Redshift Survey 
(2dFGRS) data [41] and Sloan Digital Sky Survey (SDSS) Data Release 7 which provides two 
distance measures, do.ss = rs{zd)/ Dv{0.35) = 0.1097 ± 0.0036 and do,2 = rs{zd)/ Dy (0.2) = 
0.1905 ± 0.0061 |5|, where rs{zd) is the acoustic sound horizon at the drag epoch, = 



[{1 + zyD\{z)/H{z)]^^^ . Fitting formula for z^ is defined by Eisenstein & Hu 42[. The 
xIaoi is 'Eijidi - df'){C^^){dj - df'), where di = (4.2,6^0.35), 

, 30124 -17227 , 
C-'=\ . (16) 



-17227 86977 



The fitting formula for z^ has this form: 



--- i^a6!^W-^a ^-HMn..yi. (17) 
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bi = 0.313(r]„/i2)-o-4i9 [1 + 0.607(n^/i2)0.674j ^ (^g) 
b2 = 0.238(a„/i2)°-223_ (^9) 



We also include BAO result from WiggleZ Dark Energy Survey 43|, which gives A{0.6) = 
0.452 ± 0.018. A(z) is given by 

Aiz) = (20) 
z 

The xIao, = [M^) - 0.452]V0.0182. Therefore, xIao = xIao, + x\ao,- 



D. Prior 

For the radiation, we fix = 2.469 x 10~^//i^, and the radiation energy density = 
fi^(l + 0.2271A^e//)! where ^efs is the effective number of neutrino species and is taken to 



be 3.04 [45|. We further impose the prior of = 73.8 ± 2.4 kms ^Mpc ^ from 4J]. The 
total chi-square Xtotai = X^Nia + xImb + X%ao + Xho is marginahzed over the nuisance 
parameters Qjji'^ and the reduced hubble constant h, by minimizing Xtotai with respect to 
fifo/i^ and h [45|. 



IV. OBSERVATIONAL CONSTRAINTS ON TRACKER FIELD MODELS 

The models we analyze include V{(j)) = M'^+"0-", V{(j)) = exp{Mp/(j)), V{(j)) = 
M^[exp{Mp/(f)) - 1], 1/(0) = M^exp{/3Mp/(j)), and 1/(0) = M^[exp(7Mp/0) - 1]. All a, /3 
and 7 are positive constant. The mass M is determined by requiring that the total fractional 
energy Qtotai equals to 1 in a flat universe. The initial conditions of and are arbitrarily 
chosen in the range ensuring that the scalar fleld joins the tracking solution before the last 
scattering. We solve Eq. [1] and Eq. [2] numerically in order to calculate the chi-square for 
each point in the parameter space. 

The resulting best-flt parameters for these flve models are listed in Table [B The late time 
evolution of w,f,{z) corresponding to the best-flt parameters are plotted in Fig. [21 The joint 
constraints on {flrn,C(), (^m,/^), {^m,l) sue shown in Fig. O 
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TABLE I: Fitting Results 



Model 


Best-fit parameters GoF ABIC AAIC 


Cosmological Constant 


Vim = 0.277tooi3 76.9% 


V oc (t)-"" 


Q _ n 077+0.013 

— L"-^' ' -0.013 ^ 

76.0% 6.3 2.0 

a = 0+0-070 


V cx exp(/3Mp/</>) 


^^-m 0«277 r\ n "1 

76.0% 6.3 2.0 

/? = 0+00051 


V oc [exp {-fMp/4)) - 1] 


— n 9QQ+0-015 
S'm — u.zyo_Q Q]^4 

41.0% 37.2 32.9 

^ ^ Q+0.049 a 


V oc exp(Mp/(/>) 


Vtm = 0.30ll[!;[3i^ 16.5% 57.5 57.5 


V oc [exp(Mp/0) - 1] 


Vtm = 0.307t[!:[]i^ 11.7% 65.1 65.1 



"the best-fit 7 arbitrarily approaches 



We further evaluate the GoF^ for each model (see Table [T]). The GoF gives the probability 
of obtaining data that are a worse fit to the model based on statistics, assuming that the 
model is correct. It can test the validity of a particular model. To assess the relative model 
strength, especially in comparison with the cosmological constant, we use the information 
criteria (IC). The IC are model selecting statistics encoding the tension between quality of 
fit and model complexity. They favor models that give a good fit with fewer parameters. 
The use of IC in the context of cosmological observation has been examined in 46| . In this 
paper we evaluate the Bayesian information criterion (BIG) 47|] and the Akaike information 
criterion (AIG) [48] for each model. 

The BIG is defined as BIG = —2 In Cmax + k\\iN , where Cmax is the maximum likelihood, 
which is equivalent to the minimum for gaussian errors, k is the number of parameters, 
and is the number of data points used in the fit. It comes from approximating the evidence 
ratios of models, known as the Bayes factor. A better model has a lower BIG. The AIG is 
defined as AIG = —2 In Cmax + 2/c. The AIG is derived by an approximate minimization of 
the Kullback-Leibler information entropy, which measures the difference between the true 



1 Defined as GoF=r(;//2, x^/2)/r(i^/2), where r(zy/2,x^/2) is the upper incomplete ffamma function and 
u is the degrees of freedom. 
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data distribution and the model distribution. A better model has a lower AIC The BIG 
gives stiffer penalty for extra parameters for the size of data InA^ > 2. The differences 
in BIG (ABIG) and AIG (AAIG) between each tracker field model and the cosmological 
constant are listed in Table [H 



V. CONCLUSION 



We have examined tracker field models with the potentials V{(f)) oc °, V{(j)) oc 
exp(Mp/0), V{(t>) oc exp(Mp/0)-l, ¥{(/)) oc exp(/3Mp/0) and V{(f)) oc exp(7Mp/0)-l, based 
on current observational data. It is shown that the resulting parameter constraints are strin- 
gent (see Tableland Fig.|3]). Best-fit of the two models V{(f>) oc 0"" and V{(j)) oc exp(/3Mp/0) 
are equivalent to the cosmological constant (see Tableland Fig. |2]). The best-fit of the other 
three models that do not have limits as the cosmological constant render late-time equation 
of state staying away from —1 {w^ > —0.8). The larger values of w,j) come from the intrinsic 
limits of these models. 

The poor GoF (GoF < 17%) of models V"(0) oc exp(Mp/0), 1/(0) oc exp(Mp/0) - 1 
indicates these two models are less valid. The rank of model strength is the same assessed 
either by BIG or AIG. In comparison with the cosmological constant, the three models 
V{(f)) oc exp(Mp/0), V{(j)) oc exp(Mp/0) - 1, and V{(j)) oc exp(7MpM) - 1 have ABIG > 6, 
while ABIG > 6 is consider a strong evidence against the model 46||. This shows that the 
worthiness of considering these three models, in the presence of the cosmological constant, 
is not supported by the current observational data. 

The result that both V{(j)) oc and V{(j)) oc exp{/3Mp/(j)) have the best-fit as the cos- 
mological constant suggests that other dark energy models which do not have the boundary 
w > —1 might render better fits to the data. These include the phantom models 49 1 with 

n 

w < —1, the quintom models (see [50| for a review) with w crossing —1 and the K-essence 
models (see [51| for a review). The observational constraints on these models and their 
model strength should be further studied. 

The next generation cosmological probes are expected to constrain w about ten times 
better 52|]. More stringent constraints on individual models are also expected to be obtained 



in the future (see 



53|-56| for the case of quintessence models). In the ongoing pursuit of 



revealing the nature of dark energy, the reconstruction of the general features and the model 
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based approach should be comphmentary to each other. While testing the cosmological 
constant by examining if i;; = — 1 and if w has dynamical behavior, we should also take the 
model based approach to see if there is other model worth considering. 
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FIG. 2: The late time evolution of w^{z) corresponding to the best-fit parameters. Models with 
V{(l)) oc and V{<p) oc exp{PMp/(f)) have their best-fit acting as the cosmological constant (A). 
Models with F(^) oc exp(Mp/^), F(^) oc [exp(Mp/^) - 1] and oc [exp(7Mp/^) - 1] have their 
w^{z) of the best-fit staying away from —1. 
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FIG. 3: Joint constraints on (f2„,a), (ri^)/?)) and (^^^,7). The dark gray and the Hght gray 
regions correspond to the 68.3% and 95.4% confidence regions, respectively. 
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